Edge currents in the absence of edges 
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Abstract 

We investigate a charged two-dimensional particle in a homoge- 
neous magnetic field interacting with a periodic array of point ob- 
stacles. We show that while Landau levels remain to be infinitely 
degenerate eigenvalues, between them the system has bands of abso- 
lutely continuous spectrum and exhibits thus a transport along the 
array. We also compute the band functions and the corresponding 
probability current. 



The fact that the presence of boundaries can induce a transport in a system 
with a homogeneous magnetic field is known for long |l|, || and has numerous 
consequences in solid state physics. The mentioned pioneering papers were 
followed by tremendous number of studies in which the magnetic transport 
was analyzed numerically in various models as well as experimentally. The 
obstacles at which the particle "bounces" need not be hard walls but also 
objects with openings such as various antidot lattices; for a sample of liter- 
ature see, e.g., Refs. p, ^, the following papers by the same authors, and 
references therein. The closest to the subject of this letter is a very recent 
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paper |^ where the array of obstacles is strictly one- dimensional and there 
are no boundaries to help the transport. 

One can say therefore that at the heuristic level most aspects of the 
two-dimensional magnetic transport are understood. With this fact and the 
extensive literature in mind it is a bit surprising how little attention was paid 
during a decade and a half to a strict derivation of transport properties from 
the first principles - which is after all the raison d'etre of theoretical physics. 
The situation concerning this aspect of the problem changed recently where 
several rigorous studies appeared. It was shown, e.g., that the edge currents 
in a halfplane survive a mild disorder so that away of the Landau levels the 
spectrum remains absolutely continuous 0, and that the result extends 
to planar domains containing an open wedge 

Recall that the wall producing the edge states need not be of a potential 
type. It is known, e.g., that a step of the magnetic field or another variation 
exhibiting a translational symmetry will smear again the Landau levels into 
a continuous spectrum [l^ . Similarly to the usual edge states, this type 



of propagation allows for a classical explanation in terms of the cyclotronic 



radius changing with the magnetic field - see Ref. |]1T|, Sec. 6.5. 

The aim of the present letter is to contribute to this series of rigorous stud- 
ies with a simple solvable model in which a charged quantum particle in the 
plane exposed to a homogeneous magnetic field of intensity B perpendicular 
to the plane interacts with a periodic array of point obstacles situated at the 
X axis and described by 6 potentials. We shall show that the model exhibits 
a magnetic transport of which we can with a license say that it is a purely 
quantum effect in the sense that a quantum particle propagates despite the 
fact that its classical counterpart moves on localized circular trajectories - 
apart of a zero-measure family of the initial conditions. The delta potentials 
are certainly an idealization; a more realistic model with potential-type ob- 
stacles will be discussed in a forthcoming paper. Let us mention, however, 
that point interactions play a distinguished role being the only obstacles that 
can preserve the Landau levels in the spectrum. 

Using the Landau gauge, we can write the Hamiltonian formally as 

Ha,i = i-td. + Byf - dl + ^«(5(x-xo-j£) , (1) 

where £ > is the array spacing. Since we are interested mainly in the essence 
of the effect, we use everywhere rationalized units h = c = e = 2m = 1. 
The interaction term, in particular the formal coupling constant a, needs an 
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explanation, since the two-dimensional 5 potential is an involved object. We 
adopt here the conventional definition from Ref. [1^ which determines the 
latter by means of the boundary conditions 

Li{^,aj) + 2'KaLo{ilj,aj) = 0, j = 0, ±1, ±2, . . . (2) 

with dj := {xo+ji,0), where are the generalized boundary values 

ib(x) r 

LQ{ilj,d) := lim - — — — — - , Li{i/j,d) := lim ilj{x) — LQ{i/j,d) ln|x — a| 

[x-a|^0 In |x — a| \x-a\^Ol 

(3) 

and a is the (rescaled) coupling constant; the free (Landau) Hamiltonian 
corresponds to a = cxd. Recall that since the magnetic field amounts locally 
to a regular potential in the s-wave subspace, the non-magnetic boundary 
conditions of Ref. |]12| need not be modified - see, e.g., Ref. |]13 . 



The difference between a and a reflects the nontrivial way in which the 
two-dimensional point interaction arises in the limit of scaled potentials. 
Due to the coupling constant renormalization, a caution is required when 
interpreting spectral properties of such a Hamiltonian. On the other hand, 
the two-dimensional point interaction yields a generically correct description 
of low-energy scattering which can be tested in experiments - see, e.g., a 



fresh example in |]I4|. To understand the above remark about the purely 
quantum nature of the transport here, recall the well-known expression for 
the scattering length which shows that the obstacles have a "point" character 
if e^'^" <C £, i.e., if the point interaction is strong enough. Below we shall 
show that the transport exists here for any flnite value of a. 

Using the periodicity, we can perform the Bloch decomposition in the x 
direction writing 

H^,i = — / H^A^)de, (4) 



27r 



'|<7r 



where the flber operator Ha/{9) is of the form (|^) on the strip < x < £ 
with the boundary conditions 

9;^(£-,l/) = e^^^9:^(0+,y), z = 0,l. (5) 

The Green's function of the operator Ha/{9) is given by means of the Krein 
formula App.A], 

{Hf^Ae)-zY^{x,x) = Go{x,x';6,z) 

+ (a-^(ao; 6*, z)y^Go{x, do; 9, z)Go{do, x'\ 9, z) , (6) 
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where Gn is the free Green's function and 



^{a]6^z)\= hm \GQ{a,x]d^z) ln|a?'— a|] (7) 

\x-a\^0 \ 2tX J 

is its regularized value at the point a. The Bloch conditions @ determine 
eigenvalues and eigenfunctions of the transverse part of the free operator, 

^^m{e) = + ^) vL{x) = ^ e^(--^^^)-/^ , (8) 

where m runs through integers. Then we have 

Go(f , x'; e,z) = - "'^yi^'^i^r^ Vii^m^') , (9) 

where y<,y> is the smaller and larger value, respectively, of y, y', and m^, 
are solutions to the equation 

-u"iy)+iBy+-^ + e\ u{y) = zu{y) (10) 

such that is at — oo and is at +oo; in the denominator we have 
their Wronskian. By the argument shift we get 

n , ^ ( 2-nm + d€\ 
v?^{y) = u[y^ j (11) 

and a similar relation for f^, where M,f are the corresponding oscillator 
solutions. Of course, we have W{u^,v'l^) = W{u,v). The functions u,v 
express in terms of the confluent hypergeometric functions [|1^, Chap. 13]: 

v{y) = e-^y'/^u(^^,\;By^^ (12) 

away from zero, and u is obtained by analytical continuation in the y"^ vari- 
able; together we have 



I \{y) = V^e-^y'/' 



r(^) Ti^) 



(13) 
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From here and Ref. [13, Chap. 6] we compute the Wronskian; in combination 
with (ID we get 

2(^/2B)-{3/2) /B-z\ , ,^ 

G„(^,.-9,.,.-l^r(^le— . 



^ " (^^< + 5£ J ^ (^^> + 5£ ) ^ • 



oo 

X 

m=— oo 



As expected the function has singularities which are independent of Q and 
coincide with the Landau levels, i.e., Zn = B{2n + 1), n = 0, 1, 2, . . .. Let us 
observe first that each Zn remains to be infinitely degenerate eigenvalue of the 
"full" fiber operator Ha/{0). To this end, one has to adapt the argument of 
Refs. [0, to the set of functions w^^sin (^) e-■^l'"l'/^ A; = 0, 1, . . ., with 
w := X + iy which vanish at the points of the array so the conditions @ are 
satisfied for them automatically. 

On the other hand, Ha/{6) has also eigenvalues away of Zn which we 
denote as e„(6') = el^''^\9). In view of (|^) they are given by the implicit 
equation 

a = ^{ao;e,e) (15) 
and the corresponding eigenf unctions are 

i;^^^'\x;9) = Goix,do;9,eni9)). (16) 

In order to evaluate them, we have to assess the convergence of the series in 
([T^). Using the asymptotic behavior 

{l]iy) = e^^^^^^'/' {tVBv) ^ (1 + 0{\yn) (17) 

for y =Foo, we find that the product 

27(771 + 9e\ f 2nm + 9i 



" U/< TTn — ]'"\y> + 



Bi J V B£ 
is for y ^ y' governed by the exponential term. 



exp {f {yl - yl) + {9- ^) (y> - ,<)} {\m\-' + 0{\mn) 

(18) 
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as |m| oo, while for y = y' we have 

1 



m|~^ + 0{\m\ 



so the series (|T^) is not absolutely convergent. Summing now the contribu- 
tions from ±m we see that in the limit x' ^ x it diverges at the same rate 
as the Taylor series of — (l/27r) ln( does for ( 0+. Hence we get 

^(-'^'^)= E I^^TF^-T^^^^OM^y + ^^Jj, (19) 

where C := The expression is independent of x, because the regularized 
resolvent does not change if the array is shifted in the x direction. We can 
write it by means of the first hypergeometric function alone, since 



(20) 



where ^ := VB {y + ^^^) . 

Spectral bands of our model are given by the ranges of the functions 
€„(■). Solutions of the condition (plSf ) do not cross the Landau levels, because 
^(ao; 9, ■) is increasing in the intervals (— oo, B) and {B{2n — 1), B{2n + 1)) 
and diverges at the endpoints; this is a general feature [|l^. The spectrum 
will be continuous away of z„ if the latter are nowhere constant. In view 
of the spectral condition (|15D one has to check that C,{x;6,z) is nowhere 
constant as a function of 6. Notice that each term in (|T9|) is real-analytic 
for real z and the series has a convergent majorant independent of 6; hence 
^(x; ■, z) is real-analytic as well and one has to check that it is non-constant 
in the whole Brillouin zone [— 7r/£, 7r/£). 

Suppose that the opposite is true. Then the Fourier coefficients 

/TT/e 
e(x; e, z) e'^^^ d9 (21) 

should vanish for any non-zero integer k. Since the summand in ([T9|) be- 
haves as (9(|m|~^) as \m\ oo, we may interchange the summation and 
integration. A simple change of variables then gives 



c, = - r(2C) hm / {uv) (y + — ] e''' , (22) 

-7r(2M+l) V -D 



M-»oo 
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so 

AM 



Fy{k):= Fy{^) e'"'' = , (23) 

J — oo 

where Fy{'d) := (uv) (y + -|^). The same reasoning apphes to any finitely 
periodic extension of ^(x; 6, z), hence (p3|) is vahd for each non-zero rational 
k. However, the function decays (9(|'f?|~"'^) and the integral makes sense only 
as the principal value. We shall use the above mentioned asymptotic behavior 
which implies, in particular, 

Fyi^) = 7^ + fyi^) , (24) 

where fy{'d) = 0{\{}\''^) uniformly in ?/ G [0,£]. Thus 

Fyik) = -^Koik) + fyik), (25) 



see |jT9|, 3.754.2]. Since fy G L^, the second term at the r.h.s. is continuous 



w.r.t. k and the same is then true for Fy] this means that the relation (p3[ ) 
is valid for any nonzero k. Furthermore, fy is bounded and Kq diverges 
logarithmically at A; = 0, hence J^j^ Fy{{}) e^^^ d-t} can be bounded by an 
integrable function independent of A^. Then 



/oo pN poo 

dk(l){k) lim / Fy{d)e'^^d§= / Fy{'d)^{§)d§ 
-oo N-*^J-N J-oo 

(26) 

holds for any G 5(]R), i.e., Fy{k) is the Fourier transform of Fy{d) in the 
sense of tempered distributions. Since this is a one-to-one correspondence 



20| , Thm.IX.2], we arrive at the absurd conclusion that Fy = 0. We get thus 
the following result: 

Theorem. For any real a the spectrum of H^/ consists of the Landau levels 
-B(2n+1), n = 0, 1, 2, . . ., and absolutely continuous spectral bands situated 
between adjacent Landau levels and below B. 

Let us remark that during the final stage of the work we learned about a 
similar result for a chain of point scatterers in a three-dimensional space 



with a homogeneous magnetic field [21 . Due to the higher dimensionality. 



the spectrum is purely a.c. in that case and has at most finitely many gaps. 

The above theorem says a little about the character of the transport. To 
get a better idea we solve the spectral condition (|15D numerically for several 
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INSERT FIG. 1 (spect2.eps) 

Figure 1: The eigenvalues of Ha^i{9) in the units of B. (Top) n = 20, B = 4 
(full line), B = Q (full dotted line) and B = 8 (dotted line). (Bottom) n = 1, 
B = 10 (full line) and B = 15 (dotted line). The thick lines represent the 
Landau levels. 

INSERT FIG. 2 (tokn20.eps) 

Figure 2: Probability current for n = 20, B = 4, a = 0.5 and two different 
values of 9 corresponding to the extremal points of e2o(-)- ^ = (top), 9 = 2.2 
(bottom), see the arrows in Fig. 1. The star marks the point perturbation 
position. 

values of the parameters. The results are plotted in Fig. 1 for the second and 
21st spectral band. We see that the bands move downwards with decreasing 
a and their profile becomes more complicated with the band index n; a higher 
B tends to smear the structure. 

The Bloch functions (|T^) are in general complex-valued and yield thus 

a nontrivial probability current, jn{x] 6^) = 2 Im ^^/'i"'^^(V — iA)tpi^'^^^ {x; 9). 

The current pattern changes with 9 oscillating between a symmetric "two- 
way" picture and the situations where one direction clearly prevails, these 
extremal behaviors occurring at the extrema of the corresponding band func- 
tion. This is illustrated in Fig. 2. In addition, while the pattern has pre- 
dominantly "laminar" character, in some parts current vortices may form, 
mainly in low spectral bands as it is illustrated in Fig. 3. Similar effects have 
also been observed in the numerical analysis of related models mentioned in 
the introduction. 

To sum up the above discussion, we have analyzed the behavior of a 
quantum particle in the plane exposed to a homogeneous magnetic field and 
interacting with a periodic array of point perturbations. We have shown 
that while the Landau levels survive, the spectrum develops an absolutely 
continuous part, i.e. a sequence of spectral bands. Depending on the quasi- 

INSERT FIG. 2 (toknlb.eps) 

Figure 3: Probability current for n=l,i? = 20,a = — 1 and ^ = 0. In the 
bottom graph the inset shows a vortex between the point perturbations. 
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momentum, the particle is transported along the array with zero or nonzero 
mean longitudinal momentum, and the probability current pattern may ex- 
hibit vortices in some regions. 
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Figure captions 

Figure 1. The eigenvalues of Ha,i{0) in the units of B. (Top) n = 20, 

5 = 4 (full line), B = 6 (full dotted line) and 5 = 8 (dotted line). 
(Bottom) n = 1, B = 10 (full line) and B = 15 (dotted line). The 
thick lines represent the Landau levels. 

Figure 2. Probability current for n = 20, B = 4, a = 0.5 and two different 
values of 6 corresponding to the extremal points of e2o(-)- ^ = (top), 

6 = 2.2 (bottom), see the arrows in Fig. 1. The star marks the point 
perturbation position. 

Figure 3. Probability current for n=l,-B = 20,a = — 1 and 6 = 0. In the 
bottom graph the inset shows a vortex between the point perturbations. 
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